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Abstract. The probability of a quantum particle being detected in a given solid angle is determined 
by the 5-matrix. The explanation of this fact in time dependent scattering theory is often United to the 
quantum flux, since the quantum flux integrated against a (detector-) surface and over a time interval 
can be viewed as the probability that the particle crosses this surface within the given time interval. 
Regarding many particle scattering, however, this argument is no longer valid, as each particle arrives 
at the detector at its own random time. While various treatments of this problem can be envisaged, 
C " 3 ■ here we present a straightforward Bohmian analysis of many particle potential scattering from which 

the ^-matrix probability emerges in the limit of large distances. 

Q\ • 1. Introduction 

In a scattering experiment the central quantity is the cross section, whose derivation is based on the 
. probabihty that particles are detected in certain solid angles. This probability is computed from the 

S'-matrix and is given by 

(1) J d'fci... J 
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where ^/'out(fc) is the Fourier transform of the N particle scattering state's outgoing asymptote and the 
CN ' CsfcS are the cones (in momentum space) subtending the solid angles in question. 

> ; . ^ 

I It is a common argument (see e.g. [o]) that the ^'-matrix formalism can be based on the quantum flux 

OO . given by 

fS| ■ According to ordinary statistical mechanics \ j^{x,t) ■ d(T\dt gives the probability that a particle crosses 

I the surface element dor during time dt. In [B] the connection between quantum flux and ([T|) was posed 

as a theorem, the so called flux-across-surfaces theorem. For one particle potential scattering it asserts 
that (see [Tl[3J[Tni[lI] for various versions) 



lim / f{x,t)-dcTdt^ lim / / \f{x,t)-d(T\dt= / Vout(fc) 

H ' i?.^ocJ J R^coJ J J 

' _RS flS Ce 



2 

r3 



d-'k. 



Here R is the "large" radius of a ball centered around the scattering center and E is a solid angle. RY, 
stands for the piece of the ball's spherical surface that subtends S and dcr for the positively oriented 
infinitesimal surface element. Asymptotically, the flux j"^ points outwards and in that case the flux 
integrated against the (detector) surface and over all times is the probability that the particle crosses 
that (detector) surface at some time (see e.g. [THMH])- When more than one particle is scattered, the role 
of the quantum flux is less clear [13] ■ Indeed, the crossing times are random, i.e. if more than one particle 
is scattered, one has a different crossing time for each particle and in such a setting the single-time many 
particle flux is of little relevance. Put in simple terms, what one wishes to compute is the joint probability 



(2) P (particle one crosses -RSi, particle two crosses R^2, ■•) , 
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which in the Umit i? — ?> oo should become ((T]) . 

In Bohmian mechanics the joint probabihty ^ is readily defined, where P is the measure P''' on the 
particles' initial positions with density ji/'P- By virtue of the continuity of Bohmian trajectories, first 
crossing times through boundaries of regions in space are well defined. Let tfx'^ be the first exit time of 



the lih particle from the ball xi < R then XJ{xo,t, 
time (see Figure [T|). We shall show that 
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) S M'^ is the position of the Ith particle at this 
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Figure 1. Sketch of the scattering situation for iV = 2. 

The idea for proving ^ is rather simple: Far away from the scattering center, the particles' trajectories 
become more or less straight lines directed along the asymptotic velocity 

vLixo) ■■= lim , 

which is IV'outp-distributecQ. 

The asymptotic straight line motion of the Bohmian trajectories arises because for large times the wave 
function is close to a local plane wave 

(4) e-^V«(iO-'^e^*V^out(f) , 

where this approximation needs, however, to hold in a sense stronger than L^. Establishing this is part of 
the technical work that forces us to consider only noninteracting - but nevertheless entangled - particles 
that are scattered off a fixed potential, like e.g. EPR pairs. Moreover, since the Bohmian velocity field, 
just as the flux, involves derivatives of the wave function, we shall need also stronger requirements on the 
potential than what one is accustomed to in S'-Matrix theory. 

A realistic description of the scattering process, where detection takes place at finite but random times 
has to take into account that the detection of one particle leads to a "collapse of the entangled wave 
function" which may influence the detection statistics of the remaining particles. Hence we shall prove 
^ with appropriately defined trajectories that take into account the effect of a particle's detection on 
the wave function. In this respect we note, that in Bohmian mechanics the wave function collapse does 
not happen physically but arises merely as an effective description. It is quite intriguing to see how this 
effective collapse can be handled very naturally by the so called conditional wave function, the wave 
function that arises if one or more particles' coordinates are known |121I14] . In this respect it is also 



^We use natural units where p - 
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interesting to investigate the effect of measuring devices on the wave function evolution, for example in 
an EPR experiment (see e.g. [20] )■ We have reduced tedious but straightforward technicalities as much 
as possible to better concentrate on the probabilistic reasoning, which is somewhat new in the context of 
quantum mechanics. 



1.1. Bohmian mechanics. Bohmian mechanics [^ [TTllT^IHlfTB] is a theory of particles in motion that 
is experimentally equivalent to quantum mechanics whenever the latter makes unambiguous predictions 
[12] . The state of N particles is described by their (normalized) wave function ij;{x,t), where x = 
(xi, xn) e M'^^, and by their actual configuration (positions) X ~ {Xi, X^) e M'^^. The 
wave function evolves according to the Schrodinger equation 

(5) ih^^i,{x,t) ^ H^{x,t) 

and governs the motion of the particle by (^ = 1, . , 



(6) -Xi'{x^,t)^vf(x^'{x^,t),t 




i^iX'^ixo, t),t) 



Here X'^(a;o,0) = Xq is the particles' configuration at time t — 0, mi is the mass of the Ith particle and 
V; is the gradient with respect to a;;. In ([S]) iJ is the usual Schrodinger Hamiltonian 



(7) H^-J2:^^i + Vix) = Ho + V{x) 



2mi 



1=1 

with the real valued potential V. From now on we shall use natural units h = mi = 1. 

According to Born's law, the positions of particles each having wave function ip are always |'0p-distributed 
(see [11] for a precise assertion). The distribution is equivariant: If X^(-,0) is |'0(-, 0)p-distributcd then 
X*(-,i) is |V'(-,t)p-distributed. In terms of the flow map : M^^ ^ M^^ of ®, i.e. 



[x^{xo,h))=X^{xo,h 



this means that 



(8) YH;t2) ^p^(-,ti) „ (^$f^^^^) ^ . 

For a wide class of sufficiently regular potentials and initial wave functions P'^'-almost sure global existence 
of Bohmian mechanics was proved in ^ and [22. . In particular, our setting below falls into the scope of 
Corollary 3.2 in 3 resp. Corollary 4 in [22) . 



1.2. Asymptotic behavior of Bohmian trajectories in scattering situations. In scattering situ- 
ations the Bohmian trajectories become straight lines at large distances from the scattering center. For 
non-interacting particles, 

N N N . 

V{x) = Y,Viixi) resp. H = Y,Hi = Y.(--Ai + Vi{xi)\ , 
1=1 1=1 1=1 ^ ^ 

this assertion is a direct extension (cf. [T^]) of the corresponding result for one particle [18., and will be 
used in the proof of ([3]). It holds for a large class of wave functions (see e.g. [THldH] in the one particle 
case) which we denote by G- The set Q is invariant under the time evolution e~'^' and dense in the set 
(2)/=i '^a.cXHi) of "pure" scattering wave functions. 

Theorem 1. Let V be sufficiently smooth and fast decaying and let zero he neither a resonance nor an 
eigenvalue of Hi (I — I, . . . , N ). Let tp & G with = 1. Then: 

(i) The Bohmian trajectories X^{xo,t) exist uniquely and globally in time for ¥^ -almost all initial 
configurations Xq G M'^^ . 
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(ii) For Y"^ -almost all Bohmian trajectories the asymptotic velocity vf^{xQ) := lim i^o-t) g^jgfg 

and the distribution of vf^ under has the density |V'out(')P- Moreover for all e > there exist 
T > and C > such that 



(9) 



'({ 



i3N 



(^X'^{xo,t),t)-vt{xo] 



< ct- 



yt>T\-) > 1 



}) 



Theorem [T] is based on the asymptotic structure of the wave function ip at large distances from the 
scattering center. We note that large distances imply large times. Asymptotically the wave function 
achieves local plane wave structure, meaning that dU holds in a pointwise sense and, moreover, that a 
corresponding statement holds for gradients. A crucial ingredient in proving these pointwise statements 
is the expansion in generalized eigenfunctions and the application of stationary phase methods |19) . 

2. Exit statistics 

In a realistic scattering experiments when a particle crosses or hits a detector surface some measurement 
takes place. In A^-particle scattering each particle is detected at a different time. The measurement of 
one particle produces a collapse of the entangled wave function and one needs to address the question 
whether that collapse has an effect on the detection statistics of the as yet not detected particles. To 
answer this question we shall first consider the mathematically idealized situation, where the detection 
event is simply modeled by the crossing of a particle's trajectory through a surface without invoking a 
collapse of the wave function. 



2.1. Exit statistics without collapse. The first exit time t^^{xo) of the trajectory {X^ {xQ,t), t > 0} 
from an open set A C 



1>3N 



(10) 



t^^{xo) := inf ^^t > \ Xq G A and X'''{xa,t) ^ a| 



We are interested in the sets Bi^ji := {x e 



\xi\ — xi < R}. The configuration space trajectory 



{X^{xo, t), t >0} leaves Bi^n when the Ith. particle's trajectory leaves the open ball Bn = {x e 



X < 



R}. Since X'^{xo,t) is continuous Xf {xo.tci''^ (xq)) G OBr, the sphere with radius R. Let C dBi 
be a solid angle in and R'Ei the corresponding subset of dBi^. The first exit of the Ith particle is in 
if and only if Xf {xo,tc-i'" (xo)) G RY,i. With the above we can formulate 

Theorem 2. Under the conditions of Theorem 1: 

jini^P'^da^o eK^^ I xf(a;o,tS-"(a;)) e i?!], , V? = 1, 2, . . . , 7v}) 

= jd^ki 



(11) 



d^k 



N 



where C^^ denotes the cone in R'^ subtended by the angle Sfe. 

Proof. By Theorem [1] we only need to show that 

(12) lim P'^ffa^o e K'"^ I xf (a^o, iS'" (a^)) G R^i , V/ = 1,2, 

7?.— >oo \ L 



For this we employ the asymptotic straightness of the Bohmian trajectories that holds for large times. 
Thinking of T as a large time we split the trajectory into two pieces: 



X^{xo,t)-vt(xo)t 



for < i < T , 



(13) f sup \X^'{xo,t)\+vt{xo)T 

< I Q<t<T 

X^'{xo,t)~X^'ixo,T)~vt{xo){t-T) +\X^'{xo,T)\+vtixo)T ioi- t > T . 
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Observe that global existence of the Bohmian trajectories implies that for given £ > and T there exists 
C such that 



(14) 



P'^({a;oeM3^ I sup \X^' {xo,t)\ < C]) > I ~ \ . 

^ 0<t<T ' o 



Next, by Theorem 1 (ii) for K large enough 



(15) 



Bk 



By integrating the velocities in we obtain the following estimate for the positions: There is C > 
and T > such that 



(16) 



9^' ({a^o e M^^ I |X^(aJo,t) - X^(a;o,T) - vt{x)(t - T)\ < CVt, Vt > t}) > 1 - | 



for all T > T. Observing the splitting (fT3|) we may thus conclude that there exists a such that 



(17) 



We set 



p3N\ 



\X^'ixo,t)-vt{xo)t\ <C,il + Vt), Vt>0}) >l-e. 



Gc 



t,3N\ 



\X^{xo, t) - vt{xo)t\ < Gil + Vi) , Vt > o} . 



and call the trajectories on this set asymptotically straight. Note that we have just shown that typical 

trajectories are asymptotically straight. Since the 
distribution of vf^ is absolutely continuous we may 
assume without loss of generality that ^ and 
that the cones Cs^ are open. Let vf^{xo) e Csi x 
• ■•X Cejv- Then the straight trajectory t)^;(a;o)t 
of the Ith particle crosses the surface i?S; at time 
Tfl :— R/v"^ i{xo). Thus the crossing time of the 
straight trajectory grows linearly with the distance R. 
The same is true for the distance between the point 
where i{xo)t crosses RT,i and the boundary of the 
surface i?S;. However, since the difference between 
the corresponding asymptotically straight trajectory 




Xf{xo, t) and the straight trajectory i{xo)t grows 
only sublinearly in time (~ 's/i), this difference evalu- 
ated at Tr also grows sublinearly with R. Hence, if R 
is big enough, at the distance between the point of 
crossing of the straight trajectory and the boundary 
of the surface is larger than the distance between 
Xf{xo,Tii) and i{xq)Tr. Since the trajectory is 
continuous in t this implies that X^[xQ^t) crosses 
dBji first in if and only if Axq) lies in Cs, 



FiGURE 2. Real (asymptotically 
straight) and straight trajectory. 



(see Figure [2]). 



□ 



2.2. Exit statistics with collapse. The actual detection of a particle changes the wave function; it 
collapses, so that the new spatial support of the wave function is in accordance with the measured 
position. In standard quantum mechanics such a collapse is introduced when the particle is measured 
at a given fixed time t. In the situation at hand however the detection time is random and the collapse 
must be invoked at that random time. In Bohmian mechanics the collapse of the wave function does 
not actually happen. The collapsed wave function is replaced by the well defined concept of conditional 
wave function which is based on the existence of trajectories and can be used to describe the dynamics 
of subsystems pTl[T2l[T4] . For example, - given the first particle's trajectory Xf{xQ,t) - the conditional 
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wave function for the remaining — 1 particles is 

V (X2,...,XN,t) .= 11 ^ -|p. 

The conditional wave function is random since the first particle's trajectory Xf{xo, t) is random. It 

yields the conditional probability F'^ [{xt{xo,t), . . . ,X%{xo,t)) G • | Xf (a;o,i) = x^^ = P'A™'' (■■*)(.): 
a straightforward calculation shows that 

P'^ ({xt{xo,t), X%{xo,t)) e d^xi . . . d^XN I Xt{xo,t) = 

|V;(a;i,a;2,...,a;jv,0p 3 . 2 ,3 ,3 

= — ij-T s||2 " xi . . .d XN =■■ |'0xi,t(a;2, . . . ,a;jv)| a xi...d xn 

\\^[Xi, tjll 

= pr-H;t) ^d^x,...d^XN) . 

In measurement-like situations the conditional wave function becomes the so called effective wave function 
which coincides with the collapsed wave function of standard quantum mechanics [Tll|T2l[T4] . 

Now let 

(18) T'^ixo) min{tf,i.«(a;o), t^^ H^o) ■ ■ • ,if/'"(^o)} 

be the time of first detection. Without loss of generality we assume that particle one is detected first, 
T^{xo) — tex''^(xo). Then its measured position is Xf{xo,T^{xo)) and the remaining particles' condi- 
tional wave function at time T^{xo) is 

* ""-^ ||*(xr(x.,r«(..)).....r«(x.))|| - 

To ease the discussion we assume that after its detection the particle is stuck at its measured position. 
Other more complicated scenarios where the detected particle is processed within the detector do not 
add further insight into the question at hand: eventually they would lead to the same technicalities as 
presented here. Thus the conditional wave function's time evolution after time T^(a?o) is given by 

(20) ^^"'ifa;, X.. t)^ ^-MH..+-+H.)(t-T-(.„)) H^ti^o, T^ixp)), X2,...,xj,, T^jxp)) 

' \\^p{xt{xo,T^ixo)),;--.,;T^ixo))\\ ' 

This time evolved conditional wave function then defines the evolution of the remaining (N — 1) particles 
after time T^{xo): 

d 



(21) -Xf {xo,t)^vf {Xt {xo,t),...,X% {xo,t),t), l^2,...,N,t>T''ixo) 



with "initial" conditions xf°"\xo,T^{xo)) = Xf {xo,T^{xo)). 
On the set of initial configurations Xq such that T^(xq) = tc^ '^ {xq) one is thus led to the stopped process 

X't'{xo,t) fort<T«(a;o) 



(22) Y ixo,t): >^(^x^(^^^^j.B.(^^^))^xr^'\xo,t),...,Xt;°''\xo,t)) iovt>T"{xo) 

In the following for ease of notation we denote the measure F''' conditioned on T-^{xo) — t^^'"(a;o) again 
by F'^. It is not a priori clear that the stopped process {Y^'{xo,t))t>o exists globally and uniquely for 
F'^'-almost every Xq. That this is the case is the new technical ingredient in this section. 

Theorem 3. Under the conditions and notations of Theorem\^the following holds. 

(i) For every R > the stopped process (Y^(xo,t))t>o exists globally and uniquely for -almost 
every Xq. 

(ii) Let R > 0, T > and C > 0. We define the set Vhtc of initial values Xq, which fulfill 

C 



(23) yT^ixo)>t>T: (Y''ixo,t),t) - vUxo) 
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and 



(24) 



(25) 



< 



C 



Then for all e > there exist > 0, > and > such that for all T > , C > Cg > 
and R> > 

P'^iVRTc) >l-e. 



We remark that this theorem is analog to Theorem [T] The only difference is that the error term between 
the conditional wave function and the asymptotic local plane wave (cf. (0])) and thus the error term 
between the "real" velocity vf and the asymptotic velocity v"^ ^ is now dominated by the first detection 
time T^{xo). Thus, once global existence and uniqueness of the stopped process, i.e. assertion (i), is 
established, the proof of assertion (ii) is completely analog to the proof of assertion (ii) of Theorem [1] 

(cf. m)- 

Proof of assertion (i): 

The global existence and uniqueness of Bohmian trajectories was proven in [3] under the condition that 
the wave function ip belongs to C°°{H) = Q^^pj I?(iJ"), where 'D{H") denotes the domain of iJ". Hence 
we need to show that belongs to C°°{H) where H := X]t=2 Omitting the normalization factor 

in ((H]), this amounts to proving that V(^i, O G 'D{H"') for every Xi e and n e N. Set i^„(Xi, •) := 
H"^P{Xi,-). Since C°°(-ff) and thus / \\{Hiipn){-, z)\\^d^^^''^h ^ \\Hiipn\\^ < \\H"+^i;\\^ < oo, 

R3<"-1) 

this implies \\{Hiipn){-, z)\\ < oo for almost every (with respect to Lesbesgue measure) 2; e and thus 

(Pn{-, z) £ V{Hi) ^ W^{m.^) for a. e. 2 G R^CA^-i) . 



Here H^^(R^) is the second Sobolev space. Thus we can apply an instance of the Gagliardo-Nirenberg 
inequality [TS1[T7], namely 



h^llL°=(R3) < c\\d'^u\ 



L2(R3) 



with II Z?^ 



'W||l2(R3) 



max ||Z3"u||/^2(R3) and C > independent of u e T4^^(R'^), to get 

\a\—m 



ll'^n(-,2:)|lL~(R3) < C||i:|2^</?„(-,z)||22(R3)||</?n(-,2)|li2(R3) 

Then, using Holder in the second to last step, 

\\MXi,-)\\Um--^y)= I Wn{X,,z)\^d'^''-'h< 



for a. e. z £ 



j3(JV-l) 



\Vni-,z)\\l^ 



^3(Ar-l) 



R3(«-l) R3(iV-l) 

^)lli(K3)||^„(-, ^)|li(R3)d3(^-^)z 



IPxi<^n(-,^)lll2(js.3) 



CO {■ rf3(JV-l) 

H^lA , Z)\\ ^2(Jj3-)« z 



<C \\D '^n|||_2(]g3iV-)||¥'n|||,2(]g3iV) 

for every Xi e Moreover, since \\H^pr^\\ < ||iZ"+V|| < 00, i.e. tp,, e V{H) = V{Hq) = W^{M?^), the 
last term is finite, i.e. we have just shown that indeed iy9„(Xi, •) e L^(R'^(^~^^) for arbitrary Xi € R^ 
and n e N. This proves that ■0™"'' belongs to C°° [h] 



The second step consists in proving that the dynamics of the stopped process for times t > T^{xq) is 
well defined, i.e. that for P'^-almost all initial configurations x^ the random vector 

{Xt{x^, T\x^)), • • • , X%{x,, T^ixo))) 
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belongs to the set of initial configurations Zq, for which 
(26) dT 

has unique global solutions. Let Q"^ be the set of initial conditions, for which the flow map of 
exists globally. In IB] is was shown that P'^(5'^) = 1. Set 

We wish to show that P'^'(G^) = 1 . By the definition of conditional probability 



, , B3x[0,oo) 

(27) 



'^{{Xti;t),...,X%{;t))eg^-^- I (Xf(,t),T^) = (a;i,t))P^^,^^^^^_^^^(d3:,,dO 

K3x[0,oo) 

where P^ ^ ^ (d^^idt) denotes the image measure of the random vector (X]''(-, T^), T^). Recall 

(-^1 i'tT j,T ) 

the flow map of (O . Denoting by $q j the inverse of $f g and by X/ the random variables Xi{x) — xi 

(I = 1, . . . , N), we have by equivariance that 

¥^{{Xt[;t),. . . ,X%{-,t)) eg*-^'^ \ {Xt{;T^),T^) = [x,,t)) 

=P^(-'*)((X2,...,XAr)eg'^^i- I {Xti;T^),T^)o^t.t^ix,,t)) 



for all i e M. Introducing this into (|27)) gives 
(29) 

P'A(g«)^y'p^(-.*)((X2,...,X^) e g^-^* I (Xf(,T«),T«)o$„^_, = (a;i,i))Pf^.(.^^H)^^H)('^'2;idt). 

R3x[0,oo) 

Next we observe that 

P'^(--*)((X2,...,X^) e^"^--' I (Xf(.,T«),T«)o$^,=(a,i,t)) 
(30) / , X 

= v^-i-' (g^-i-' I r^^ o = t]^i 

where the second equality is a direct consequence of P'^^'i ' (fj'/'^^i *) = 1. Introducing this in (I^S)) we get 
the desired almost sure global existence P'^(G^) = 1 . To see the first equality in ([50]) we use 

and compute: Let D C x and split x = ixi,z). Then 

P'^(-'*)(Zeg'^-i'%(Xi,T«o$^Jei?) = y" |i^(a;,i)|2d'^a; 

(||V'(a;i,-,i)||2 y |^,^,t(^)pd'('^'')^)rf'a;i 

{zee'*-i.'|(a!i,T«o*;f Jxi,z))e_D} 

\\yj{x,,;t)\\^P^-^^* {g^-^^\T^ o ^tAx,,-) e {s I (a;i,s) e i^D^^^i 



(|iV(a.i,.,t)ir y P^'-'(t;^-'MT«o$;f^,(a.i,.) = .)Pj;;-;.^(^^_.)(ds))d3:ri 

{s\(xi,s)eD} 



D 



i.e. by the definition of the conditional probability 

pV-C-,*) (^z e g'^-i'* I (Xf(-,T^),T^) o = (a;i,s)) = P^-i'' (g^'-i^* | o = 
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for all s, t e M. Setting s — t gives the first equality in (15(11) . 

3. Conclusion 

It is of course desirable to extend the results to iV-particle scattering, i.e. to the case of interacting 
particles. However, our analysis here relies heavily on earlier works, where it was shown that the scattered 
wave function attains local plane wave structure Q with a small pointwise error. This was achieved 
using an expansion into generalized eigenfunctions, for which regularity properties had been established 
in |10ll21j . To apply this method to interacting particles we would need similar regularity properties for 
the generalized eigenfunctions also in this case. These are however not known and presumably much 
more difficult to get. 
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